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Abstract. We recall how parton distributions are constructed in a statistical physical picture 
of the nucleon. The chiral properties of QCD lead to strong relations between quarks and 
antiquarks distributions and the importance of the Pauli exclusion principle is also emphasized. 
A global next-to-leading order QCD analysis of unpolarized and polarized deep-inelastic 
scattering data allows to determine a small number of free parameters. Some predictions are 
compared to recent experimental results and we discuss the prospects of this physical framework. 



1. Introduction 

Deep-inelastic scattering (DIS) of leptons on hadrons has been extensively studied, over the 
last twenty years or so, both theoretically and experimentally, to extract the polarized parton 
distributions (PPD) of the nucleon. As it is well known, the unpolarized light quarks (n, d) 
distributions are fairly well determined. Moreover, the data exhibit a clear evidence for a flavor- 
asymmetric light sea, i.e. d > u, which can be understood in terms of the Pauli exclusion 
principle, based on the fact that the proton contains two u quarks and only one d quark. Larger 
uncertainties still persist for the gluon (G) and the heavy quarks (s,c) distributions. From the 
more restricted amount of data on polarized structure functions, the corresponding polarized 
gluon and s quark distributions (AG, As) are badly constrained and we just begin to uncover a 
flavor asymmetry, for the corresponding polarized light sea, namely Am ^ Ad. Whereas the signs 
of the polarized light quarks distributions are essentially well established, An > and Ad < 0, 
this is not the case for An and Ad. In this talk, essentially based on Refs.^|2J|3], we will recall 
the construction of a complete set of polarized parton (all flavor quarks, antiquarks and gluon) 
distributions and, in particular, we will try to clarify this last point on the polarized light sea. 
Our main motivation is to use the statistical approach to build up : q, L , Aqi, q, L , Aqi, G and AG, 
by means of a very small number of free parameters. A flavor separation for the unpolarized and 
polarized light sea is automatically achieved in a way dictated by our approach. Our predictions 
are compared with some new unpolarized and polarized DIS measurements obtained at DESY, 
SLAC and Jefferson Lab., as well as hadronic cross sections, from FNAL and RHIC-BNL, and 
they turned out to be very satisfactory. 
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2. Construction of the PPD in the statistical approach 

In the statistical approach the nucleoli is viewed as a gas of massless partons (quarks, antiquaries, 
gluons) in equilibrium at a given temperature in a finite size volume. The light quarks q = u,d 
of helicity h = ±, at the input energy scale Ql = 4GeV 2 , are given by the sum of two terms pQ, 
a quasi Fermi-Dirac function and a helicity independent diffractive contribution, common to all 
light quarks 

h 2 _ AX 0g xb Ax 1 

xq (x, Q ) = exp[( ^ _ - i + exp(x/ - } - 1 • (1) 

Here Xq is a constant, which plays the role of the thermodynamical potential of the quark q h and 
x is the universal temperature, which is the same for all partons. Clearly one has qi = qf + q~ , 
for the unpolarized quark distributions and Aqi = qf — q[ for the polarized ones and similarly 
for antiquarks and gluons. We recall that from the chiral structure of QCD, we have two 
important properties, allowing to relate quark and antiquark distributions and to restrict the 
gluon distribution: 

- The potential of a quark q h of helicity h is opposite to the potential of the corresponding 
antiquark q~ h of helicity -h 

^Oq = ~^-0q 

The potential of the gluon G is zero 



X 0q - ~ X 0a ■ ( 2 ) 



X 0G = . (3) 

Therefore similarly to Eq. l[T|). we have for the light antiquarks 

h 2 A{X^)- l x 2b Ax 1 

% ~ q (X ' Qo) = expKx + V)/^ 1 + ^ x l~ x ) + 1 ' ^ 

Here we take 2b for the power of x and not b as for quarks, an assumption which was discussed 
and partly justified in Ref. pQ. 

Concerning the unpolarized gluon distribution, we use a quasi Bose-Einstein function, with zero 
potential 

xG(x,Ql)= f GX _" G . (5) 
exp(x/x) — 1 

This choice is consistent with the idea that hadrons, in the DIS regime, are black body cavities 
for the color fields. It is also reasonable to assume that for very small x, xG(x, Qq) has the same 
behavior as the diffractive contribution of the quark and antiquark distributions in Eqs. (1) and 
(4), so we will take be = 1 + b. We also need to specify the polarized gluon distribution and we 
take 

xAG(x,Q 2 ) = , (6) 

assuming a zero polarized gluon distribution at the input energy scale Qq . 

For the strange quarks and antiquarks s and s, given our poor knowledge on their unpolarized 

and polarized distributions, we take 

xs{x, Ql) = xs(x, Ql) = - [xu(x, Ql) + xd(x, Q§)] , (7) 



and 



xAs(x, Ql) = xAs(x, Ql) = \ [xAd(x, Q 2 ) - xAu(x, Ql)} . (8) 



This particular choice gives rise to a large negative As(x,Qq). Both unpolarized and polarized 
distributions for the heavy quarks c, b, t, are set to zero at Qq = 4GeV 2 and for larger Q 2 , they 
are generated by the Q 2 evolution. 

With the above assumptions, we note that the heavy quarks do not introduce any free 
parameters, likewise the gluons, since the normalization constant Aq in Eq. (5) is determined 
from the momentum sum rule. Among the parameters introduced so far in Eqs. (Q) and (4), 
A and A are fixed by the two conditions u — u = 2, d — d = 1. Clearly these valence quark 
conditions are independent of b and A, since the diffractive contribution cancels out. Therefore 
the light quarks require only eight free parameters, the four potentials X^ u , Xq u , X^ d , X^ d , one 
universal temperature x and three more parameters b, b and A. 

From well established features of the u and d quark distributions extracted from DIS data, 
we anticipate some simple relations between the potentials: 

- u(x) dominates over d(x), therefore one expects X^ u + X^ u > X^ d + X^ d 

- Au(x) > 0, therefore X^ u > X^ u 

- Ad(x) < 0, therefore X~ d > X+ d . 

So Xq u should be the largest thermodynamical potential and X^ d the smallest one. In fact, 
as we will see below, we have the following ordering 

Ku > X od ~ X 0u > K d ■ (9) 

This ordering leads immediately to some important consequences for quarks and antiquarks. 
First, the fact that X^ d ~ Xq u , indicated in Eq. (9), leads to 

u-(x,Q 2 )<d-(x,Q 2 ) , (10) 

which implies from our procedure to construct antiquark from quark distributions, 

u + {x,Ql) >d + (x,Q 2 ) . (11) 



These two important approximate relations were already obtained in Ref. fj]. by observing in 
the data, the similarity in shape of the isovector structure functions Ixgf n \x) and F^f n \x), 

at the initial energy scale, as illustrated in Fig. 1. For Ixgf n \x) the black circles are obtained 
by combining SLAC jl] and JLab data. The white circles, which extend down to the very 
low x region, include the recent deuteron data from COMPASS 6 combined with the proton 
data from SMC [7j, at the measured Q 2 values of these two experiments. The agreement with 
the curve of the statistical model is improved in this later case. The + helicity components 
disappear in the difference 2xg^f n \x) — n \x). Since this difference is mainly non-zero for 
0.01 < x < 0.3, it is due to the contributions of u~ and d~ (see Ref. U). 

Second, the ordering in Eq. (9) implies the following properties for antiquarks, namely: 

i) d(x) > u(x), the flavor symmetry breaking which also follows from the Pauli exclusion 
principle, as recalled above. This was already confirmed by the violation of the Gottfried sum 
rule [HI EH] - 

ii) Au(x) > and Ad(x) < 0, which have not been established yet, given the lack of 
precision of the polarized semi-inclusive DIS data, as we will see below. One expects an 
accurate determination of these distributions from the measurement of helicity asymmetries 
for weak boson production in polarized pp collisions at RHIC-BNL ^H] , which will allow 
this flavor separation. By performing a next-to-leading order (NLO) QCD evolution of 
these parton distributions, we were able to obtain in Ref. pQ, a good description of a large 
set of very precise data on the following unpolarized and polarized DIS structure functions 



F% ,d,n (x, Q 2 ) : xF^ N (xj Q 2 ) and g^' ,n (^ ? Q 2 ), in a broad range of x and Q 2 , in correspondance 
with the eight free parameters : 

X+ u = 0.46128, X Qu = 0.29766, X~ d = 0.30174, X+ d = 0.22775 , (12) 

x = 0.09907, 6 = 0.40962, b = -0.25347, A = 0.08318 , (13) 
and three additional parameters, which are fixed by normalization conditions 

A = 1.74938, A = 1.90801, A G = 14.27535 , (14) 

as explained above. Note that the numerical values of the four potentials are in agreement with 
the ordering in Eq. (9), as expected, and all the free parameters in Eqs. (12, fT3|) have been 
determined rather precisely, with an error of the order of one percent. 
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Figure 1. The isovector structure functions 
Ixgf n \x) (solid line from our statistical parton 

distributions) and F^ (x) (dashed line from our 
statistical distributions). Data are from SLAC [I], 
JLab 0, COMPASS @, SMC 7 and NMC [H] 
(Taken from Ref. [3]). 



Figure 2. F^x^Q 2 ) as function of Q 2 for 
fixed x, c(x) = 0.6(i x — 0.4), i x = 1 — ► x = 
0.32, rebinned data HI, ZEUS, E665, NMC, 
BCDMS. (Presentation of data, courtesy of 
R. Voss) (Taken from Ref. pQ) . 



3. Experimental tests for unpolarized and polarized DIS 

We first consider (j,p and ep DIS for which several experiments have yielded a large number of 
data points on the structure function F% (x, Q 2 ). We have compared our calculations with fixed 
target measurements NMC, BCDMS and E665, which cover a rather limited kinematic region 
in Q 2 and also with the data at HERA from the HI and ZEUS Collaborations. These last data 
cover a very large Q 2 range, up to Q 2 = 10 4 GeV 2 or so and probe the very low x region which 
is dominated by the rising behavior of the universal diffractive term. We compare our results 
with the data on Fig. 2. We also have a very good description of the neutron structure function 
F^ix, Q 2 ) data, as well as for the differential inclusive neutrino and antineutrino cross sections 
from the high statistics vN DIS, CCFR and NuTeV experiments as shown on Figs. 3 and 4 . 
As expected, for fixed x, the y dependence is rather flat for neutrino and has the characteristic 
(1 — y) 2 behavior for antineutrino. 
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Figure 3. Differential cross section v(v)N 
for E v = 85GeV and 0.045 < x < 0.225, 
as a function of y. Data are from CCFR 
(white circles) and NuTeV experiments 
[El IE] (black circles) (Taken from Ref. 0). 



Figure 4. Differential cross section v(v)N 
for E v = 85GeV and 0.275 < x < 0.65, 
function of y. Data are from CCFR 
11 (white circles) and NuTeV experiments 
[Ellin] (black circles) (Taken from Ref. 0). 



Since our approach is based on the direct construction of the quark and antiquark 
distributions of a given helicity and qf , from the previous results we immediately obtained 
Aqi and Aq~i for each flavor, which enter in the definition of the polarized structure functions 
g^ ,d,n (x , Q 2 ) . In Fig. 5 we show a data compilation of g\ )d)n {x,Q 2 ) from different current 



experiments on proton, deuterium and helium targets, evolved at a fixed value Q 2 = 5GeV 2 . The 
x dependence is in fair agreement with our results and we predict, in the small x region, a fast 
rising behavior for g\ and a fast decreasing behavior for g™, due to the antiquark contributions. 
This cannot be tested so far, due to the lack of precise data. The very high x region also 
remains poorly known and the typical behaviour of the Fermi-Dirac functions, falling down 
exponentially above the thermodynamical potential, which shows up in Fig. 1, complies well 
with the fast change in the slope of g\{x) at high x, as shown in Fig. 6. 
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Figure 5. gi' d,n (x,Q 2 ) versus x for different 
Q 2 values, from E155, E154, E143, SMC, 
HERMES. The curves correspond to our 
model predictions at Q 2 = 5GeV 2 (Taken 
from Ref. HI). 



Figure 6. g\{x,Q 2 ) 3>S EL function of x 
at fixed Q 2 = 3GeV 2 from the statistical 
approach. Experimental data from SLAC 
E143 dH (Taken from Ref. 0). 



Let us now turn to the helicity distributions, whose best determinations are shown in Figs. 7 
and 8. The ratio Au(x)/u{x) is predicted to have a rather fast increase in the x range 
(Xq u — x,Xq u + x) and a smoother behaviour above, while Ad{x)/d{x), which is negative, 
has a fast decrease in the x range (X^~ d — x, X^ d + x) and a smooth one above. This is exactly 
the trends displayed in Fig. 7 and our predictions are in perfect agreement with the available 
accurate high x data. We note the behavior near x = 1, another typical property of the statistical 
approach, also at variance with predictions of the current literature. The fact that Au(x) is more 
concentrated in the higher x region than Ad(x), accounts for the change of sign of g™(x), which 
becomes positive for x > 0.5, as first observed at Jefferson Lab |16| . 

Concerning the antiquark helicity distributions, it is clear that in the very low x region 



(x < 10 3 ) the ratio d(x)/u(x) is ~ 1, because the diffractive contribution dominates. Since we 
predict Au(x) > and Ad(x) < 0, it is also interesting to remark that Eq. (11) implies 



Au(x) — Ad(x) ~ d(x) — u(x) > 



(15) 



so the flavor asymmetry of the light antiquark distributions is almost the same for the 
corresponding helicity distributions. Similarly, Eq. (10) implies 



Au(x) — Ad(x) ~ u(x) — d(x) > . 



(16) 



We also compare in Fig. 8 our predictions with an attempt from Hermes to isolate the different 
quark and antiquark helicity distributions. The poor quality of the data does not allow to 
conclude on the signs of Au(x) and Ad(x), which will have to wait for a higher precision 
measurement of helicity asymmetries in production at RHIC-BNL |181 119j . 
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Figure 7. Ratios (Au + Au)/(u + u) 
and (Ad + Ad)/(d + d) as a function of x. 
Data from Hermes for Q 2 = 2.5GeV 2 [E] 
and a JLab experiment |16j . The curves 
are predictions from the statistical approach 
(Taken from Ref. 0). 



Figure 8. Quarks and antiquarks polarized 
parton distributions function of x for 

Q 2 = 2.5GeV 2 . Data from Hermes [13. 
The curves are predictions from the statistical 
approach (Taken from Ref. [B]). 



4. Single jet and ir inclusive productions 

A precise determination of parton distributions allows us to use them as input information to 
predict strong interaction processes, for additional tests of pertubative QCD and also for the 



search of new physics. Here we shall test our statistical parton distributions for the description 
of two inclusive reactions, single jet and tt° productions. The cross section for the production 
of a single jet of rapidity y and transverse momentum pr, in a pp collision is given by 



E^ = Et-V-/ 1 ^^^ 
dp 6 *-ri 1 + dij 7r J XQ 2x a - xj 

% 3 



XTG y 

fi(xa,Q 2 )fj(xb,Q 2 )^r(s,i, u) + {i <r-yj) 
at 



(17) 



where xt = 2pt/^/s, xq = XTe y /(2 — XTe~ y ), Xf, = x a xxe~ y / (2x a — XTe y ) and ^fs is the center 
of mass energy of the collision. In the above sum, i,j stand for initial gluon-gluon, quark-gluon 
and quark-quark scatterings, d&ij/dt are the corresponding partonic cross sections and Q 2 is the 
scaling variable. 

In Fig. 9 our results are compared with the data from CDF and DO experiments |2U1 121j . 
Our prediction agrees very well with the data up to the highest Et (or pr) value and this is 
remarkable given the fact that the experimental results are falling off over more than six orders 
of magnitude, leaving no room for new physics. For completeness, we also show in Fig. 10 the 
DO data, for several rapidity bins, using a (Data-Theory) /Theory presentation. 
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Figure 9. Cross section for single jet 
production in pp at = 1.8TeV as a 

function of Et- Data are from CDF j^H] and 
DO [U experiments (Taken from Ref. 0). 



Figure 10. Comparison between the 
statistical model and the DO H2] single jet 
cross sections in pp at yfs = 1.8TeV, as a 
function of Et and rapidity r/ (Taken from 
Ref. 0). 



Next we consider the cross section for the inclusive production of a tt° of rapidity y and 
transverse momentum px, in a pp collision, which has the following expression 

E n d 3 a/dpl = J dx a dx b f a / p (x a ,Q 2 ) x 

abc 

fb/ P (x b , Q ) 1 ' da/dt{ab - cX) , (18) 

irz c 

where the sum is over all the contributing partonic channels ab — ► cAT and dcr /di is the associated 
partonic cross section. In these calculations f a /p,fb/p are our parton distributions and D^j c 
is the pion fragmentation function. Our calculations are done up to the NLO corrections, 
using the numerical code INCNLL of Ref. [22j and for two different choices of fragmentation 
functions namely, BKK of Ref. [23] and KKP of Ref. [2^], and we have checked that they 
give similar numerical results. We have compared our predictions to two different data sets at 
y/s = 200GeV from PHENIX and STAR at RHIC-BNL. The results are shown in Figs. 11 and 
12 and the agreement is good, both in central rapidity region (PHENIX) and in the forward 
region (STAR). This energy is high enough to expect NLO QCD calculations to be valid in a 
large rapidity region, which is not the case for lower energies |28| . 
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Figure 11. Inclusive 7r° production in pp 
reaction at yfs = 200GeV function of 
Pt, scale fx = px- Data from PHENIX 
|26j . Solid curve fragmentation functions 
from KKP |25j . dashed curve from BKP 
H3] (Taken from Ref. 0). 



Figure 12. Inclusive 7r° production in pp 
reaction at yfs = 200GeV as a function of 
E w . Data from STAR j^Tj. Solid curve 
fragmentation functions from KKP |25j . 
dashed curve from BKP [23| (Taken from 
Ref. 0). 



5. Final remarks 

We have shown that this simple approach of the statistical parton distributions, which involves 
only eight free parameters, provides a good description of recent data on unpolarized and 
polarized DIS and on some hadronic processes. Our distributions have special features, which 
remain to be tested by new precise experimental data in the future, as well as its natural 
extension to the transverse degree of freedom |29j . 
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